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We describe an experiment in which we have used a cold 
damping feedback mechanism to reduce the thermal noise of a 
mirror around its mechanical resonance frequency. The moni- 
toring of the brownian motion of the mirror allows to apply an 
additional viscous force without any thermal fluctuations as- 
sociated. This scheme has been experimentally implemented 
with the radiation pressure of an intensity-modulated laser 
beam. Large noise reductions, up to 30 dB, have been ob- 
tained. We have also checked the mechanical response of the 
cold damped mirror, and monitored its transient evolution 
between the cooled regime and the room temperature equi- 
librium. A simple theoretical model allows to fully explain 
the experimental results. A possible application to the active 
cooling of the violin modes in a gravitational-wave interfer- 
ometer is discussed. 

PACS : 05.40.Jc, 04.80.Nn, 42.50.Lc 



I. INTRODUCTION 

Thermal noise is a major limitation to many very sen- 
sitive optical measurements such as interferometric grav- 
itational wave detection jl],^|. The suspended mirrors 
of interferometric detectors have many degrees of free- 
dom: pendulum modes of the suspension system (with 
resonance frequencies below the analysis band), violin 
modes of the suspension wires (with some resonances 
within the detection band), and mirrors' internal acoustic 
modes (with resonance frequencies above the band) ||. 
The thermal motion of each degree of freedom results 
from its excitation by a thermal random force, called the 
Langevin force. The fluctuation-dissipation theorem re- 
lates the spectrum of the Langevin force to the temper- 
ature and to the mechanical damping of the mode. Up 
to now, only passive methods such as the increase of me- 
chanical quality factors [|| or cryogenic methods to lower 
the temperature ^ have been considered to lower the 
thermal noise and increase the sensitivity. 

It has been recently proposed to use a feedback loop to 
reduce the thermal noise of a mirror in a high-finesse op- 
tical cavity 0. A Fabry-Perot cavity is actually a very 



sensitive displacement sensor and a feedback mech- 
anism can be used to control the brownian motion of 
a mirror, for example via the radiation pressure of an 
intensity-modulated laser beam This approach al- 
lows to reduce the thermal noise both at the mechanical 
resonance frequencies of the mirror and at low frequency. 
At resonance, this cooling corresponds to a cold damping 
mechanism, the radiation pressure of the light applying a 
viscous force to the mirror without any additional ther- 
mal noise. 

The purpose of this paper is to present a more de- 
tailed theoretical analysis of the experiment described in 
||, focussing on the resonant case, and to present new 
experimental results on the mechanical response of the 
cold damped mirror and on the transient behavior of the 
system. 

Section |l| is dedicated to a brief reminder of the 
fluctuation-dissipation theorem and its consequences on 
the thermal noise spectrum. Section [II presents a de- 



tailed theoretical analysis of the cold damping process 
used in section IV to lower the effective temperature of 
the fundamental acoustic mode of a mirror. In section |v| 
we take into account the background thermal noise due 
to the other acoustic modes of the mirror to explain the 
discrepancy between experimental results and the simple 
monomode theory of section [II. 

We present in section VI the full mechanical character- 
ization of the cold damped mirror. We have studied the 
mechanical response of the mirror to an external force, 
when the cooling is applied. The results show that the 
mechanical susceptibility is changed in agreement with a 
cold damping mechanism: the additional radiation pres- 
sure corresponds to a viscous force and there is no ad- 
ditional noise associated with this force, except the elec- 
tronic noise of the feedback loop. 

Finally, we present the transient evolution of the ther- 
mal noise when the cold damping is switched on or off 
(section VII). The results show that the relaxation time 



towards the cooled regime can be much shorter than the 
relaxation time towards the thermal equilibrium. These 
transient characteristics may be used to perform a cyclic 
cooling of the thermal noise associated with the violin 
modes of the mirrors in a gravitational- wave interferom- 
eter. 
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II. THERMAL NOISE AND 
FLUCTUATION-DISSIPATION THEOREM 

We present here some well known results on thermal 
noise and on the fluctuation-dissipation theorem, which 
will give some physical insight to the efficiency of the 
feedback scheme used to lower the temperature. 

Each degree of freedom is equivalent to a harmonic 
oscillator, whose thermal motion is the consequence of 
the \ksT thermal energy and appears as the mechanical 
response of the oscillator to a Langevin force Ft describ- 
ing the coupling to a thermal bath. In the framework of 
linear response theory |J, the resulting motion can be 
described by its Fourier transform Sx [ft] at frequency ft 
which is proportional to Ft [ft] , 



Sx [ft] = \ [ft] Ft [ft] 



(1) 



where \ [ft] is the mechanical susceptibility of the mode, 

1 
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(2) 



with a mass M, a resonance frequency ftjvfj and a loss 
angle </>(ft). In the following, we consider only the case 
of a viscous damping, e.g.: 



(3) 
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where Y is the mechanical damping of the mode. It is well 
known that such a loss angle may not be appropriate to 
describe some degrees of freedom, such as internal acous- 
tic modes of the mirrors [ flO|| : <j> IS usually considered as 
independent of frequency over a large frequency range. 
However, as we are only interested in the spectral analy- 
sis of the motion at frequencies close to a mechanical res- 
onance, for low loss oscillators {(f) "C 1), both descriptions 
are equivalent. The former, leading to a simple physical 
interpretation, will be used throughout the paper. 

The spectrum Sf t of the Langevin force is related by 
the fluctuation-dissipation theorem J9| to the mechanical 
susceptibility, 



S Ft [ft] = Im 



1 



xM 



(4) 



According to Eqs. (Q) to (^), the noise spectrum S% [ft] 
of the displacement 5x at temperature T has a lorentzian 
shape, 



sj[n] =| x [fi]| 2 s Ft [v] 

2Yk B T 



(Q 2 M - n 2 y + r 2 ft 2 



(5) 



Performing the integration of (Q) over frequency yields 
the thermal variance Axt, which obeys the equipartition 
theorem, 



-MQ 2 M Ax 2 



(6) 



If we see the thermal motion as a response to the 
Langevin force Ft, it should be clear that increasing the 
damping would reduce the thermal variance. However the 
damping Y cancels from that expression, because it ap- 
pears both in the mechanical susceptibility and in the 
spectrum of the Langevin force. Decreasing the damping 
however concentrates the thermal noise around the me- 
chanical resonance frequency, and therefore lowers both 
the low frequency and high frequency thermal noises, 
which is the reason for the ongoing works on low loss 
materials for the substrates of the mirrors and the sus- 
pension wires [Q. 

The total noise spectrum, taking into account all the 
mechanical modes, is the sum of such lorentzian compo- 
nents. Near one of the mechanical resonances, the ther- 
mal noise is mainly ruled by the resonant component, 
with a background noise related to all other modes and 
approximately flat in frequency. 



III. COLD DAMPING THEORY 

The cooling mechanism is based on a measurement of 
the thermal noise of the mirror and on a feedback loop 
which applies a properly adjusted force on the mirror. 
We will see in section [V a practical implementation of 
such a mechanism. Let us only assume here that we are 
able to measure the position Sx [ft] of the mirror, and to 
apply a viscous feedback force proportional to the 
speed —iQSx [ft] of the mirror, 



F fb [ft] = -MTg x (-iQSx [ft]) , 



(7) 



where g is a dimensionless parameter related to the gain 
of the feedback loop. In this section we assume that the 
mirror can be described as a single harmonic oscillator. 
The resulting motion is given by 



Sx [ft] =x[ft] (F T M+F fb [ft]), 



(8) 



where \ [ft] is the mechanical susceptibility (0) with a 
loss angle of the form (|j|), and Ft is the Langevin force 
whose spectrum is given by the fluctuation-dissipation 
theorem (^|). The motion can be written as 



Sx [Vt]=Xfb [ft] [ft] 



(9) 



where Xfb is an effective mechanical susceptibility given 
by 



Xfb [ft] 



M(ft| f - ft 2 - i(l + g)Yfl)' 



(10) 



The additional viscous force (|7|) has obviously increased 
the damping, but, unlike the case of passive damping, 
the Langevin force is not modified and still verifies the 
fluctuation-dissipation theorem (||): it is only related 
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to the natural damping T of the oscillator. We will 
see shortly that this so called cold damping mechanism 
jll],[l2| allows to reduce the effective temperature of the 
system. 

Equations (^) and ( |l0| ) show that the mirror now re- 
sponds to the Langevin force with the effective mechan- 
ical susceptibility Xfb [^]- As the spectrum of Ft is flat 
against frequency, the thermal noise spectrum in pres- 
ence of feedback S£ [tt] still has a lorentzian shape, of 
width (1 + g)T, 



s[ b [n] = 



2Tk B T 



1 



M 



{fl 2 



n 2 y + (i + g ) 2 r 2 n 2 



(ii) 



As the random force has not been increased, the variance 
Ax 2 of the motion is reduced. Performing the integration 
of (Thy) over the frequency yields 



-Mn 2 M Ax 2 



k B T 
2(1 + 5)' 



(12) 



The cold damped oscillator of damping T and at thermo- 
dynamical temperature T is therefore equivalent to an 
oscillator with an effective damping 



r /6 = (i + s )r, 

at an effective temperature 

T fb = T/(l + g). 



(13) 



(14) 



This can also be deduced from the spectrum of the 
Langevin force which can be written as 



S Ft = 2MTk B T = 2MT fb k B T fbl 



(15) 



where both Tff, and Tf b have already been defined. The 
increase of the effective damping therefore decreases the 
effective temperature of the oscillator. The cold damped 
mirror is still at thermodynamical equilibrium but, unlike 
the case of passive damping, the temperature is reduced. 
We can also define the amplitude noise reduction 



R[Q] 



SIM 

Sl b [Q] 



n 2 



Q 2 



ft 2 



n 2 - im 



(16) 



Off resonance, R [f2] goes to unity: as the mechanical sus- 
ceptibility gets real, the change of its imaginary part does 
not contribute anymore to a modification of the spectrum 
S[ b [f2] . The amplitude noise reduction is maximum at 
resonance where it is equal to 



R T = R{VLm] = 1 



(17) 



We have shown in this section that measuring the po- 
sition of an oscillator allows the realization of a feedback 
loop to apply an additional viscous force to the oscillator 
without any additional thermal fluctuations. Such a cold 
damping scheme allows to set the oscillator at thermal 
equilibrium at a lower temperature. This is the case onl 
for a single harmonic oscillator. We will see in section 
that there is no longer a thermodynamical equilibrium 
when one takes into account all the acoustic modes. 



IV. OBSERVATION OF THE COLD DAMPING 

In this section, we present the experimental implemen- 
tation of the feedback scheme presented in the previous 
section. The oscillator is the fundamental acoustic mode 
of a plano-convex resonator. We present the experimen- 
tal results on the cooling at resonance, and the depen- 
dence of the effective damping and temperature with the 
gain of the feedback loop. 



A. Presentation of the oscillator and motion sensor 

The oscillator used throughout the paper is the funda- 
mental acoustic mode of a mirror coated on the plane side 
of a plano-convex resonator made of fused silica, whose 
thermal noise observation has been reported in J/j. The 
mirror has a thickness (at the center) of 1.5 mm, a diam- 
eter of 14 mm, and a curvature radius of the convex side 
of 100 mm. With the value of sound velocity in silica, 
the value of the thickness yields a fundamental resonance 
a bit below 2 MHz. 
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FIG. 1. Experimental setup. The brownian motion of the 
plano-convex mirror is measured by a high-finesse cavity. A 
frequency (F. stab.) and intensity (I. stab.) stabilized laser 
beam is sent into the cavity and the phase of the reflected field 
is measured by homodyne detection. The radiation pressure 
of an intensity-modulated beam is used to excite the motion 
of the mirror, either to characterize its mechanical response 
or to cool it with a feedback loop. 

The mirror has been coated at the Institut de Physique 
Nucleaire de Lyon and is used as the end mirror of a 
single-port high-finesse optical cavity (figure |l|) . The cou- 
pling mirror is a Newport high-finesse SuperMirror, held 
at 1 mm of the rear mirror. The light entering the cavity 
is provided by a titane-sapphire laser working at 810 nm 
and frequency locked to an optical resonance of the high- 
finesse cavity. The light beam is also intensity-stabilized 
and spatially filtered by a mode cleaner. 

The phase of the field reflected by the cavity is very 
sensitive to changes in the cavity length. It is monitored 
by a homodyne detection. The signal is superimposed 
to the quantum phase noise of the reflected beam. With 
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our parameters (finesse of the cavity T ~ 37 000, wave- 
length of the light A ~ 810 nm and incident power 
Pin — 100 fjW), the quantum limited sensitivity is on 
the order of 2.7 x 10~ 19 m/^/Hz 0. The quantum noise 
is therefore negligible and the phase of the reflected field 
reproduces the brownian motion of the mirror. 




FIG. 2. 



1858.4 1858.6 1858.8 

Frequency (kHz) 

Phase noise spectra of the field reflected by the 
cavity for a frequency span of 500 Hz around the fundamen- 
tal resonance frequency of the mirror. Curve a reflects the 
brownian motion of the mirror at room temperature. Curve b 
is the theoretical thermal noise deduced from the mechanical 
respons e of t he mirror to the optical excitation (curve c; see 
section 



IVB) and from the fluctuation-dissipation theorem. 



Curve (a) of figure g shows the phase noise spectrum 
of the reflected field for frequencies around the funda- 
mental resonance frequency of the mirror. The noise is 
calibrated in displacements of the mirror and reflects the 
brownian motion of the mirror which is peaked around 
the resonance frequency. This experimental result allows 
to deduce the mechanical resonance frequency of the fun- 
damental acoustic mode of the resonator: 



2tt x 1858.6 kHz. 



(18) 



The mechanical susceptibility has a lorentzian shape with 
a width 



r ~ 2tt x 43 Hz. 



(19) 



In other words, the mechanical quality factor Q = CIm /T 
has a value of 44 000 which seems to be limited by clamp- 
ing losses. The effective mass of the mode is deduced 
from the equipartition theorem, 



M ~ 230 mg. 



(20) 



B. Implementation of the feedback loop 

To apply a controlled force on the mirror, we use an 
auxiliary laser beam reflected from the back on the mirror 
(see figure |l|) . This beam is derived from the same laser 
source as the one used for motion sensing. The beam 
is intensity-modulated by an acousto-optic modulator so 
that a modulated radiation pressure is applied to the 
mirror: 



F rad (t) = 2hk I (t) , 



(21) 



where k — 2tt/X is the wave vector of the laser and I (t) 
is the intensity modulation, expressed as the number of 
photons per second. F rac i therefore appears as the mo- 
mentum transferred to the mirror during the reflection of 
a single photon, multiplied by the number of reflections 
per second. 

Driving the acousto-optic modulator with a monochro- 
matic excitation then allows to characterize the dynamics 
of the mirror at that frequency. We scan the excitation 
frequency for 50 values spanned around the resonance fre- 
quency in order to characterize the mechanical response 
of the mirror to an external force, as this has been done in 
. The mechanical response shown in curve (^c) displays 
a resonance centered on the same frequency as the ther- 
mal noise, and with the same width. Moreover, the ob- 
served thermal noise spectrum is in excellent agreement 
with the spectrum Sj [Ci] deduced from this mechanical 
response and the fluctuation-dissipation theorem (curve 
||b) |Q. This confirms that the observed thermal peak is 
related to the brownian motion of the fundamental acous- 
tic mode of the plano-convex resonator at room temper- 
ature. 

Note that the optical excitation is very efficient: with 
an approximately 100 mW completely modulated auxil- 
iary beam, one gets a displacement at resonance 30 dB 
larger than the one associated with the thermal noise for 
a 1 Hz resolution bandwidth of the spectrum analyzer. 

This setup can be used to implement the feedback loop 
discussed in section [II . Indeed, the signal 5q ou t [CI] given 
by the homodyne detection is proportional to the dis- 
placement of the mirror: 



Sq ou t [Ci] — a[Ci]Sx [Ci] 



(22) 



where the factor a [Ci] depends on the analysis frequency 
CI through the filtering by the cavity. In the following, 
we will focus on the thermal noise spectrum on a 1 kHz 
band around the fundamental resonance frequency. This 
band is small compared to the bandwidth of the cavity 
(on the order of 2 MHz), so that a [CI] can be consid- 
ered constant. Moreover, this band is small compared 
to the mechanical resonance frequency CIm, so that CI 
can also be considered constant and the derivation neces- 
sary to obtain the force (0) from the displacement Sx [CI] 
can be practically obtained by a dephasing of the homo- 
dyne signal. The feedback loop is therefore simply made 
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of an amplifier with adjustable gain and phase, in or- 
der to drive the acousto-optic modulator with a signal in 
quadrature with the displacement Sx [ft] . 

The signal is also filtered by an electronic bandpass fil- 
ter inserted in the loop in order to avoid the saturation 
of the feedback. Indeed, the homodyne signal also com- 
prises resonances associated with higher order modes of 
the plano-convex resonator (at higher frequencies) and 
resonances of the input mirror (mainly at lower frequen- 
cies). Moreover, even for the fundamental acoustic mode 
considered here, the spectrum of the Langevin force is 
flat against frequency and may lead to the saturation of 
the loop. 

As we are only interested in the spectrum around the 
mechanical resonance frequency, the filter is centered on 
this frequency. Its bandwidth is chosen large enough in 
order to limit its effect on the frequency band of 1 kHz 
around the resonance, and small enough in order to ef- 
ficiently filter the low-frequency and high-frequency tails 
of the resonance. The results presented here have been 
obtained with a 10 kHz bandwidth, which fulfills both 
conditions. The influence of the filter can therefore be 
neglected in our theoretical model. 

C. Noise reduction at resonance 

Figure || presents the thermal noise spectra obtained 
without feedback, and with feedback for increasing val- 
ues of the gain g of the loop m. The thermal noise at 
the mechanical resonance frequency is strongly reduced. 
Meanwhile, as predicted by our theoretical model, the 
width of the resonance is strongly increased. The effec- 
tive temperature of the oscillator is related to the vari- 
ance Ax 2 and thus to the area of the curve. The decrease 
of these areas thus corresponds to a cooling of the mirror. 
We will see in the following that large cooling factors can 
be obtained. 




Frequency (kHz) 



FIG. 3. Noise spectra obtained without feedback (curve a) 
and with feedback for increasing values of the gain (curves 
6). The temperature, proportional to the area of the ther- 
mal peak, is decreased. The 300 Hz shift in the resonance 
frequency relative to figure |i| is due to a room temperature 
drift. 



D. Dependance of the cooling with the gain 

In order to check our theoretical model, it is instruc- 
tive to study the efficiency of the cooling mechanism with 
respect to the gain g of the feedback loop. To measure 
the feedback gain, we detect the intensity of the auxil- 
iary beam after its reflection on the mirror. According 
to Eq. (j^), the ratio between the modulation spectrum 
of the intensity at frequency VLm and the noise spectrum 
S[ b [fljvf] of the displacement 6x leads to the gain <?, a 
multiplicative factor aside. Such a measurement takes 
into account any nonlinearity of the gain due to a possi- 
ble saturation of the acousto-optic modulator. 

Figure |] shows the relative damping Tt — T fb /T and 
the amplitude noise reduction at resonance Rt experi- 
mentally observed. These parameters are derived from 
the experimental spectra by lorentzian fits. As expected 
from Eqs. (|l^) and (|l7|), both have a linear dependence 
with the gain. The straight line in figure ^ is in excellent 
agreement with experimental data, and allows to nor- 
malize the gain <?, as this has been done in the figure. 
However, the cooling factor T/Tfb, derived from the area 
of the lorentzian fit, does not obey the linear dependence 
(|l4|). The following section is dedicated to the explana- 
tion of this behavior. 




Gain g 

FIG. 4. Variation of the relative damping Ft (squares), of 
the amplitude noise reduction Rt at resonance (circles), and 
of the cooling factor T/Tfb (diamonds), as a function of the 
feed back gain g. Solid lines are theoretical results (see section 
VA|). 



V. INFLUENCE OF THE BACKGROUND 
THERMAL NOISE 

The purpose of this section is to derive a more so- 
phisticated theory of our experiment, in order to explain 
the nonlinear dependence of the effective temperature 
with the gain of the loop. This dependence is due to the 
background thermal noise associated with other acous- 
tic modes of the mirror or with the input mirror of the 
cavity. 



5 



A. Dependence of the temperature with the gain 

We can take this background noise into account in our 
theoretical model. Suppose we note Sxb [ft] the associ- 
ated Fourier component, and Sb [ft] its noise spectrum. 
As there is no acoustic resonance around the frequency 
band of interest, Sb is supposed to be frequency indepen- 
dent. The noise dxb [ft] is uncorrelated to the Langevin 
force Ft [ft] of the fundamental mode, and the overall 
noise spectrum without feedback is the sum of the reso- 
nant and background components, 



Si [ft] = Sj<°> [ft] + S b 



(23) 



where Sx [ft] is the noise spectrum of the fundamental 



mode (Eq. |5j). 

Moreover, as the feedback does not involve the input 
mirror, and as we have already seen that the efficiency 
of the loop strongly decreases off resonance, we can state 
that the background thermal noise is unaffected by the 
loop. We then have 

Sx [ft] = x [ft] {F T [ft] + F fb [ft]) + 5x b [ft] , (24) 

where \ [^] is the mechanical susceptibility of the funda- 
mental acoustic mode and the radiation pressure Ffb is 
related to the overall displacement Sx by Eq. (Q). The 
resulting motion is 



Sx [ft] = X fb [«] F T [ft] 



ftL - ft 2 - irft r 

Sx h 



ft 2 



ft 2 -ir fb £l 



[ft] 



(25) 



where Xfb [ft] is the effective susceptibility of the funda- 
mental mode given by Eq. (0). As Sx b [ft] and F T [ft] are 
uncorrelated, the noise spectrum takes the simple form: 



Sl b [ft] = 5f(°) [ft] 



i?[ft] 



rSb, 



(26) 



where sl b ^ [ft] is the lorentzian spectrum obtained when 
neglecting the background noise (Eq. p"l| ), and R [ft] is 
the amplitude noise reduction defined in Eq. (|l6|). Com- 
bining Eqs. (16) and (p6j) leads to 



s f x b m 



stb 
i?[ft] 2 



(27) 



The feedback loop still decreases the noise spectrum by 
a factor R [ft] 2 . In other words, the effect of the feedback 
is unchanged since it is still described by the mechanical 
susceptibilities \ an d Xfb of the fundamental mode only, 
but the background noise in S*J [ft] leads to a background 
noise evenly modified in S[ b [ft] . The noise spectrum 
can also be written as 



Sl b [ft] = (l- £b . 9 (2 + 5 ))Sf(°) [Q] + S b 



(28) 



where Eb — Sb/ Sx^ [ftiw] is the ratio between the back- 
ground noise spectrum and the resonant component at 
frequency ftjv/ of the thermal noise without feedback 
(Eq. |23|). The overall noise spectrum stays lorentzian, the 
area of the lorentzian component being smaller when tak- 
ing Sxb [ft] into account, with an unchanged background 
noise. 

The width of the resonance is still equal to Tji, = 
(l+g)T, and the overall amplitude noise reduction is 
unchanged and equal to R [ft] as we have experimentally 
observed in figure ||. 

The effective temperature Tfb of the fundamental 
acoustic mode, related to the area of the resonant compo- 
nent of the spectrum by the equipartition theorem, can 
be computed using Eq. 



1 



T 



T fb 1 - e b g(2 + g) 



(29) 



The background thermal noise (e& > 0) thus decreases 
the effective temperature. 

One can experimentally obtain the ratio Eb using the 
lorentzian fit of the thermal noise spectrum observed 
without feedback (figure ^|a). The experimental value 
is Eb — 1/150. The solid curve in figure [| is a fit of 
the experimental points using Eq. (p9|). The optimum 
value is Eb — 1/110, close to the expected value. The 
theoretical curve is in excellent agreement with the ex- 
perimental points: taking the background thermal noise 
into account thus allows to explain the dependence of the 
effective temperature with the gain of the feedback loop. 




1858.6 1859 
Frequency (kHz) 

FIG. 5. Noise spectra obtained without feedback (curve a) 
and with feedback for a very large gain (curve b). One gets 
a dip in the background thermal noise (dashed line) and the 
mirror is no longer in a thermodynamical equilibrium. 

The effective temperature Tfb defined previously holds 
for the resonant mode, but as the background is un- 
changed by the feedback loop, one can say that all the 
other acoustic modes stay at room temperature T. The 
mirror is thus no longer in a thermodynamical equilib- 
rium. For sufficient large gains (g > 10 in our exper- 
iment), the feedback loop can dig a hole in the back- 
ground noise, leading to a negative effective temperature 
for the resonant mode. Figure || shows the thermal noise 
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spectrum experimentally observed for g ~ 40. This fig- 
ure also shows that very significant noise reduction can 
be obtained: the noise power is reduced at resonance by 
a factor larger than 1000. 



B. Background noise and thermal equilibrium 

We have discussed in the previous section how the 
background thermal noise alters the cooling mechanism. 
We have used a phenomenological model, where the back- 
ground thermal noise is described as an additional noise 
5xb [A] which is unaffected by the feedback loop. If we 
neglect the contribution of the coupling mirror to the 
background noise, it is possible to derive the same result 
more formally. 

One can show that the motion of the mirror as seen 
by the intracavity light can be described by an effective 
mechanical susceptibility which takes into account all the 
acoustic modes and their spatial overlap with the laser 
beams 13: 



Xeff [A] 



[n] (u n ,v ) 



(30) 



where \n [A] is the mechanical susceptibility of the acous- 
tic mode n with spatial deformation u n (r) along the axis 
of the cavity at every point r of the surface, vq (r) is the 
spatial profile of the intensity for the intracavity field 
and for the cooling beam (we assume that both beams 
have the same profile) . (/, g) here stands for the spatial 
overlap between two two-dimensional functions, 



(f,g) 



/(r) 5 (r) d 2 r. 



(31) 



The thermal motion is then described by an effective 
Langevin force, which is a sum of the Langevin forces 
associated with each acoustic mode. One can show |T|] 
that this Langevin force obeys the fluctuation-dissipation 
theorem with respect to Xef f > and the thermal noise spec- 
trum is therefore given by 



Sl[£l] = \Xeff [A] i 



! S Fr [<l] = ^Im( X e//[fi]). (32) 



Around any resonance of the mirror, most of the dynam- 
ics is ruled by the resonant mode, and Xeff [A] mostly 
displays the resonant behavior of this mode. 

We can now compute the motion of the mirror with 
the feedback. The motion is as previously 



6x [A] =Xfb [A] F T [A] 



(33) 



where Xfb [A] is the mechanical susceptibility of the cold 
damped mirror, 



Xfb [A] = - 



f 



I xeff [si]-iMvsi g ' 



(34) 



Using Eqs. (|3J) and 
trum takes the form: 



the displacement noise spec- 



Si h [A] 



Xfb [A] 



[(>■ [A] 

Xeff ["] 

VXeff [A] 



1/Xeff [A] - iMTQg 



SI [A] 



(35) 



As a first approximation, one can replace the effective 
susceptibility Xef f [A] around the fundamental resonance 
of the mirror by its resonant component, and find the 
thermal noise spectrum in presence of feedback: 



Sl b [A] 



f 



R[Sl] 



I Si [A] 



(36) 



where R [A] has already been defined (see Eq. . As 
previously, we find that the overall thermal noise spec- 
trum is reduced by a factor R [A] 2 . The corrections in- 
volved by the off resonance terms in Xef } [A] are negligi- 
ble for our experimental conditions: the existence of all 
other acoustic modes leads to the presence of a back- 
ground thermal noise in the noise spectra S"J [A] and 
S[ b [A] , but the dynamics of the cooling mechanism is 
mainly described by the resonant component of Xef f [A] . 

This model also allows to show that the mirror is no 
longer in thermodynamical equilibrium in presence of 
feedback. Indeed, this would require to be able to de- 
fine an effective temperature Tfj, obeying the fluctuation- 
dissipation theorem at every frequency A: 



S Ft [A] = - 



2k R T 



fb 



A 



lm(l/ X fb [A]) 



(37) 



Using Eqs. ( |34| ) and (p7|), one can show that this requires: 



Tfb = T/ 



1 - 



MTClg 



Im(l/ Xe // [A]) 



(38) 



This condition can be fulfilled only if Im (1/xe// [A]) de- 
pends linearly of the frequency A. This is the case if 
Xef f [A] describes the mechanical response of a single 
harmonic viscously-damped oscillator, but it is no longer 
the case for a mechanical resonator with many acoustic 
modes. 



VI. 



MECHANICAL RESPONSE OF THE COLD 
DAMPED MIRROR 



The experimental results presented in section IV have 



shown that it is possible to reduce the observed thermal 
noise at resonance by a factor (1 + g) with the feedback 
loop. However the observed quantity is the dephasing of 
the reflected field, which drives the feedback loop. This 
error signal goes to zero as the gain of the loop is in- 
creased. In other words, the fact that the thermal noise 
is reduced does not prove that the feedback corresponds 
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to a cold clamping mechanism. To check if the residual 
brownian motion of the mirror is given by Eqs. (|p) and 
© (or Eqs. II and H in presence of a background ther- 
mal noise), we can determine the cold damped mechani- 
cal susceptibility Xfb [ft] by measuring the mechanical re- 
sponse of the mirror to an external force. This, combined 
with the observation of the brownian motion, allows to 
check that the Langevin force Ft is not altered by the 
feedback loop, and that the residual brownian motion is 
given by Eq. (|). 



A. Response to an external force 

The motion of the mirror under the combined influ- 
ence of the Langevin thermal force Ft, the cold damping 
force Ffb, and the external force F ext used to measure its 
mechanical susceptibility is 

Sx [ft] = Xeff [ft] {F T [ft] + iMTflg Sx [ft] + F ext [ft]) 

= X fbp]{F T p] + F ext [Q}) 1 (39) 

where the effective mechanical susceptibility Xfb [ft] has 
been defined in Eq. |34|). In presence of feedback, the 
mirror responds the same way to the Langevin force and 
to an external force. Checking the mechanical response, 
one must observe a widening and an amplitude decrease 
as the gain g of the loop is increased. 



B. Experimental results 

The cold damped mechanical susceptibility Xfb [ft] is 
measured by the same technique already used in jf). 
We use a second auxiliary beam, which is intensity- 
modulated by an independent acousto-optic modulator 
controlled by a high-frequency oscillator. The two auxil- 
iary beams are mixed by a polarizing beamsplitter before 
their reflexion on the mirror. The intensity of the second 
beam is 50 mW, corresponding to a 30 dB displacement 
at resonance with respect to the thermal noise level. It is 
lower than the one of the cooling beam (500 mW, 50 dB) 
in order not to alter the cooling mechanism. 

For each value of the gain g of the loop, we measure 
the thermal noise spectrum, which is used to calibrate g. 
We then proceed to a characterization of the response of 
the mirror to the external force in presence of feedback. 
The response is measured for 50 values of the modulation 
frequency of the second auxiliary beam, spanned around 
the mechanical resonance. Figure |^ represents the mod- 
ulation power measured for different values of the gain 
of the loop. Note the similarity between this figure and 
figure H 

Each of the curves of figures || and |^ can be compared 
to a lorentzian fit, taking into account the presence of a 
background thermal noise. One gets the relative width Tt 
of the thermal resonance, the relative width Tm of the 



modulation resonance, the amplitude thermal noise re- 
duction Rt at resonance (defined in Eq. |l7|) , and the am- 
plitude modulation reduction Rm at resonance obtained 
from the ratio between curves (b) and (a) of figure ||. 
Figure ^ represents the dependence of these parameters 
as a function of the gain of the loop (Ft is not displayed 
for g > 10 because the width of the noise spectrum Sl b 
is then too large to be accurately measured with a 1 kHz 
acquisition bandwidth) . 




1858.5 



1859 
Frequency (kHz) 



FIG. 6. Mechanical response of the mirror to the exter- 
nal force without feedback (curve a) and with feedback for 
increasing values of the gain (curves b). Solid lines are 
lorentzian fits. 




FIG. 7. Variation of the relative dampings Ft (squares), 
Tm (stars), of the amplitude noise reductions at resonance 
Rt (triangles), Rm (circles), for the thermal noise and for 
the mechanical response, as a function of the feedb ack gain g. 
Solid lines are theoretical results (see section VI C| ). 



For g < 30, the experimental points are well aligned 
on a unity slope straight line. This confirms that the 
parameters obey the expected linear dependence. The 
susceptibility Xfb [ft] perfectly describes the mechanical 
behavior of the cold damped mirror. Moreover, the ther- 
mal noise spectrum, characterized by Tt and Rt, still 
obeys the fluctuation-dissipation theorem with the me- 
chanical susceptibility Xfb [ft]- As stated in our theoret- 
ical model, the Langevin force Ft [ft] is not altered by 
the cooling mechanism. 
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C. Influence of the electronic noise of the loop 

However, for sufficient large gains, the measured ther- 
mal noise reduction Rt is lower than expected with our 
theoretical model. This is due to the electronic noise of 
the feedback loop. We can take this noise into account in 
our model. Suppose we note 5x e [ft] the electronic noise 
of the loop. The cold damping force then takes the form 



F fb [ft] = iMTSlg (Sx [ft] + Sx e [ft]) , 



(40) 



and the resulting motion is 

Sx [ft] = Xfb [ft] (Fr [ft] + iMYVtg Sx e [ft] + F ext [ft]) . 

(41) 

The response to an external force F ext [ft] is not altered 
by the electronic noise, and it is still described by the ef- 
fective susceptibility Xfb [ft] - However, the observed noise 
without any external force does depend on the electronic 
noise 5x e [ft]. As the electronic noise is uncorrelated to 
the Langevin force, and flat in frequency on the frequency 
band of interest, the thermal noise spectrum appears as 
the sum of the thermal noise S[ h [ft] in presence of feed- 
back (see Eq. and the electronic noise spectrum S e 
injected by the feedback loop, 



s x [ft] = sl b [ft] + \Mm gX fb [ft] 1 2 s e . 



(42) 



The thermal noise spectrum therefore still has a 
lorentzian shape of width (1 + g)T, as it is proportional 
to \\fb [ft] | , but its amplitude does depend on the elec- 
tronic noise. The thermal noise at the mechanical reso- 
nance frequency is 



S x [ftAf] — 



1 



9 



(43) 



The influence of the electronic noise therefore increases 
with the gain g of the loop. The amplitude noise reduc- 
tion at resonance is now 



Rt 



1 



9 



(44) 



Ml 



The electronic noise can be experimentally measured 
from the characteristics of the amplifier used in the loop, 
and one gets S e /S x [SIm] = 2.2 10~ 4 . The solid curve in 
figure ^ shows the expected dependence of Rt with the 
gain, without any adjustable parameter. It is in excellent 
agreement with experimental data. 



VII. TRANSIENT REGIME 

We have focussed up to now on the steady-state of 
the cold damped mirror. Some questions however arise. 
What is the timescale of the transition from the natural 



regime (at thermodynamical temperature T = 300 K) to 
the cold damped regime (at effective temperature T/b)? 
And when the cooling mechanism is switched off? It is 
well known that the quality factor Q (and therefore the 
damping T) can be measured through the decay time of a 
mechanical excitation abruptly switched off. Is it still true 
for the cold damped regime? In this section, we present 
the measurements of the timescale needed for the mirror 
to reach its equilibrium state, first when the cooling is 
applied or switched off abruptly, and then in the case of 
a mechanical excitation, both for the open-loop and the 
cold damped cases. 



A. Transition to and from a damped regime 

Suppose we apply sharply the cooling mechanism at 
t = 0. For t > 0, the equation of motion is 



M [x + Tx + Q 2 M x] = F T - MFgx, 



(45) 



where Fr(t) is the temporal evolution of the Langevin 
force. The last term of Eq. (|4^) represents the cold 
damping force. One has to solve this equation with ini- 
tial conditions related to the previous thermodynamical 
equilibrium at temperature T. The effect of the cooling 
mechanism can be monitored on the variance Ax 2 (t), 
which is directly related to the effective temperature by 
the equipartition theorem. We show in Appendix A that 
the time evolution of the variance is given by 



Ax 2 {t) = 



1 + fle 



-r fh t 



(46) 



where Ax?r is the thermal variance without feedback. The 
temperature therefore decreases to its new equilibrium 
value Ax\j (1 + g) with a time constant rjb- ^ n other 
words, the mirror evolves towards its damped regime 
with a mechanical response corresponding to its effective 
susceptibility Xfb [ft] in presence of feedback. 

In the opposite case, when the cooling is abruptly 
switched off, the variance obeys (see Appendix A): 



Ax z (t) = Ax 2 T 1 - 



,-rt 



1+.9 



(47) 



The variance goes back to its thermal equilibrium value 
Axt with a time constant T _1 which is now related to 
the mechanical damping without feedback. The transient 
evolution to the cooled regime is then faster by a factor 
Tfb/T = 1 + g than the relaxation towards the thermal 
equilibrium. 

To observe these transient regimes, the evolution of the 
variance is monitored with the spectrum analyzer, cen- 
tered on the resonance frequency in the zero-span mode. 
A fast electronic switch is inserted in the loop in order to 
switch alternatively the loop on and off. We have checked 
that the transition time of the switch is negligible as com- 
pared to the mechanical time-scales. 







The resolution bandwidth of the analyzer must be large 
enough in order to reduce the temporal filtering caused 
by the spectrum analyzer, but not too large compared 
to the width of the mechanical resonance, so that the 
background noise does not contribute too much to the 
observed signal. With a bandwidth of 1 kHz, both con- 
ditions arc fulfilled, the time constant of the analyzer 
being on the order of 0.2 ms. 




Time (ms) 

FIG. 8. Time evolution of the displacement noise power (in 
arbitrary units) when the cooling is applied (r on ) and when 
it is switched off (t //). Dashed curves are theoretical fits. 

Figure || shows the evolution of the variance. When 
the cooling is applied, the variance decreases with a 
time constant r on ~ 0.5 ms, and when it is switched 
off, the variance relaxes towards its equilibrium value 
with a time constant r ff ~ 2.6 ms. The experimental 
data are well fitted by Eqs. ( f46| ) and (Ejj), with time 
constants in reasonable agreement with the expected 
value: T ~ 2ir x 43 Hz leads to a theoretical value 
r ff ~ 3.7 ms. The gain of the loop is here g = 5.2, 
so that the ratio T ff/r on = 5.2 is also in reasonable 
agreement with its expected value of 1 + g = 6.2. 



B. Decay time for a mechanical excitation 

This paragraph is dedicated to a direct measurement 
of the damping of the resonance, through a measure- 
ment of the decay time after a mechanical excitation, 
both with and without feedback. When a monochro- 
matic force F ext — Fq cos (SIm^) is abruptly applied to 
the mirror at t = 0, we show in Appendix B that the 
displacement spectrum (without feedback) becomes 



S X [Q 



M 



F 



4 \MSl 



M 



)'( 



1 



-rt/2 



(48) 



when the thermal noise is neglected as compared to the 
effect of the modulation. If the monochromatic excitation 
is switched off at t = 0, then the displacement spectrum 
is (see Appendix B): 



So: [SI 



Ml 



1 ( F 



4 \MSl 



M 



-rt 



(49) 



In presence of feedback, we have already seen that the 
cold damped mirror is equivalent to a mirror of intrinsic 
damping Tfj, at temperature X/&. As the thermal noise 
is neglected, we get similar expressions for S x [SIm] i n 
presence of feedback, replacing T by Tfi,. 

Figure ^| shows the observed spectra when the ex- 
citation is alternatively switched on and off, without 
feedback. The modulation power is measured as pre- 
viously by the spectrum analyzer in zero-span mode 
with a resolution bandwidth of 1 kHz. The mechani- 
cal excitation is applied by the second auxiliary beam 
intensity- modulated at frequency SIm (see section VI B) 
and switched on and off by the fast electronic switch. The 
experimental data are well fitted by Eqs. (48) and flls|), 
with a time constant of 2.8 ms, in excellent agreement 
with the one found in the previous paragraph. 
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Time (ms) 

FIG. 9. Time evolution of the modulation power of dis- 
placement (in arbitrary units) when the external force is 
switched on or off, without feedback. Dashed curves are the- 
oretical fits. 
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Time (ms) 

FIG. 10. Time evolution of the modulation power of dis- 
placement (same unit as in figure ^) when the external force 
is switched on or off, with a feedback gain of 5.2. Dashed 
curves are theoretical fits. 

Figure |l^ shows the same spectra for a cold damped 
mirror, with a gain g = 5.2. The two time constants 
deduced from the fits with Eqs. ( [i"8| ) and ( |49| ) are 0.5 ms 
and 0.6 ms, both of the same order of magnitude with 
the one found in the previous paragraph in the study of 
the relaxation towards the damped regime. 
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We have therefore shown that the relaxation times of 
the cold clamping mechanism are the ones related to the 
mechanical response to an external force, with or without 
feedback. 



switched off. The gravitational wave still has the same 
effect on the interferometer, whereas the thermal noise is 
reduced: one then gets a gain in the signal to noise ratio. 



VIII. CONCLUSION 

We have demonstrated that it is possible to actively re- 
duce the thermal noise of a mirror with a feedback loop. 
The cold damping mechanism described in this paper has 
allowed us to obtain large noise reductions around the 
resonance frequency of the fundamental acoustic mode 
of the mirror, up to 30 dB. A simple theoretical model, 
taking into account both the background thermal noise of 
the other acoustic modes of the mirror and the electronic 
noise of the loop, allows to explain quantitatively the ex- 
perimental results. We have also studied the mechanical 
response of the cold damped mirror, and the transient 
regime of the cold damping. 

We have in particular studied the mechanical response 
of the cooled mirror to an external force, and we have 
shown that the cooling mechanism indeed corresponds 
to a cold damping, where a viscous force is applied on 
the mirror without any additional thermal noise but the 
electronic noise of the feedback loop. 

We have also studied the transient regime of the cold 
damping mechanism. We have related the relaxation 
times when the cooling is switched on or off to the me- 
chanical decay times in response to an external force. Our 
results show that the cooled regime is much faster to es- 
tablish than the thermal equilibrium once the feedback 
is switched off. 

The cold damping might be used to perform a cyclic 
cooling of thermal peaks associated with mechanical res- 
onances of the mirror. One can indeed cool a mechanical 
resonance during short periods of time, leaving the sys- 
tem unaffected during much longer periods. One then 
gets a non stationary thermal noise since it is quickly 
reduced during the beginning of the cycle and slowly in- 
creases the rest of the time. If the period of the cycle 
is properly chosen as compared to the relaxation time of 
the resonance, the noise will stay well below its thermal 
equilibrium value. Furthermore, the duration of the cool- 
ing phase can be very short if the gain of the feedback 
is large enough: the system is therefore free most of the 
time. 

Such a cyclic cooling may be used to detect signals 
corresponding to single bursts of duration smaller than 
the period of the cycles. For example, one can envision 
to cool the violin modes in a gravitational-wave interfer- 
ometer which have sharp resonances within the detection 
band. Quality factors of these modes are very large so 
that the period of the cyclic cooling can be made much 
longer that the expected duration of gravitational waves. 
For a large feedback gain, the most probable situation is 
then to detect a gravitational wave when the cooling is 
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APPENDIX A 

In this appendix, we study the relaxation of the ther- 
mal variance Ax 2 (t) between two equilibrium states. In 
the case where the cooling is applied at t = 0, one has 
to solve equation ( f45|) with initial conditions correspond- 
ing to the thermal equilibrium at temperature T . The 
Laplace transform of (W5ft takes the form 



[ P 2 + r fbP + n 2 M ] x [ P ] = + (r + P ) x (o) + x (o) , 



M 



(50) 



where x [p] and Ft [p] are respectively the Laplace trans- 
forms of the displacement x (t) and of the Langevin force 
Fr(i). x(0) and £(0) are the initial conditions of the 
oscillator. This leads to 

x [p] = Xfb [p] (F T [p] + M (r + p) x (0) + Mx (0)) , (51) 

where Xfb [p] is t ne Laplace transform of the effective 
susceptibility, 



Xfb [p] 



1 



M[n 2 I+ p 2 + r fbP Y 



(52) 



x (t) can be obtained from Eq. (|5l| ) by inverse Laplace 
transform. It appears as the sum of the response x p (t) 
to the Langevin force (first term in Eq. glj) and the 
transient evolution xq (t) from the initial conditions 
M0),±(0)}: 



with 

x F (t) 
x Q (t) 



x (t) — x F (t) + x (t) 



Xfb (t-r)F (t) dr, 



(53) 



(54a) 



M (Tx (0) + x (0)) Xfb (t) + Mx (0) Xfb (t) , (54b) 



where the temporal evolution of Xfb (t) is obtained by 
inverse Laplace transform of (j52|): 



Xfb (t) 



sin (Cl M t) 



M 



• exp 



(55) 
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As Ft (t > 0) and the initial conditions are uncorre- 
lated, the variances of xf and xq must be summed, 



Ax 2 (t) = Ax 2 F (t) + Ax 2 , (t) . 



(56) 



The Langevin force corresponds to a white noise, there- 
fore its correlation function is 



F T (t) F T (f) = 2MTk B T S (t - t') , 
and the variance Ax 2 ? (t) is equal to 

Ax% (t) = 2MTk B T f [ Xfb (r)] 2 dr. 
Jo 



(57) 



(58) 



In the limit of a large quality factor of the resonance 
(r <c £Im), this leads to 



Ax% (t) = 



Axi, 



1 



1 



-T fb t\ 



9 



(59) 



where Ax T = k B T/Mfi 2 M is the variance at room tem- 
perature. 

At t — 0, the mirror is in thcrmodynamical equilib- 
rium at temperature T, thus we get Ax 2 (0) = A; 



Ax 2 (0) = Of^ Axf, and x (0) x (0) = 0. Combined with 
the expression of xq (t) (Eq. |54b| ) , this leads to 



Axl (i) = Ax T e~ Vlbt . 
The variance Ax 2 (t) is finally equal to 



Ax 2 (t) 



At 2 



1 



3e 



-r /b t 



(60) 



(61) 



One can derive in a similar way the result when the 
cooling is switched off at time t = 0, by taking g = 
in Eq. ( [45| ) and initial conditions corresponding to the 
cooled regime. One gets 



Ax 2 (t) 



Aa4 | 1 



.9 _ e -rt 



1 + 3 



(62) 



x(t) 



X(t-r) F ext (r) dr, 



(64) 



where X (^ — r ) is the same as in Eq. (55), except for 
the damping T. For a monochromatic force F ex t (t) = 
F cos (ilAjt + 4>)i one gets 



x(t) 



F 



MQ M T 



sin (fijiff + 



1 - e- rt / 2 



(65) 



As r ^ Om , we have a slow increase of the amplitude of a 
monochromatic spectrum. In that limit, one can compute 
the temporal evolution of the displacement spectrum at 
resonance, 



4 \MO M r 



1 - e- rt / 2 



(66) 



In the case where the monochromatic excitation is 
switched off at t — 0, one has to take into account the ini- 
tial conditions which correspond to the forced regime of 
the oscillator. The Laplace transform of the displacement 
obeys 



[ P 2 +T P + n 2 M ] x [p] = (r + p) x (0) + x (0) , (67) 



where: 



:(0) = 



F 







MQ M T 



sin (</)) . 



Ft 



x (0) = TJf cos (^) • 
Using the results of Appendix A, one finds 
x(t) =Mx(0) X (t)+Mx(0)x(t) 



(68a) 
(68b) 



F n 



■ sin (Qm t 



-rt/2 



MS] M r 

and the corresponding displacement spectrum is 



s x [n 



M 



F 



1 



4 V Mfl M T 



-rt 



(69) 



(70) 



APPENDIX B 



This appendix is dedicated to the study of the transient 
response of the mirror when an external force is switched 
on or off. Suppose we apply at t — an external force 
Fext to the mirror, without feedback. The displacement 
x(t) then obeys 

M [x + Tx + n 2 M x] = F ext + F T . (63) 

If the response to the external force is much larger than 
the thermal noise, we can neglect the Langevin force and 
the initial conditions, which are related to the brownian 
motion of the mirror. One then gets 
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